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STRESSES  IN  GUN  TURKS 
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The  stress  distribution  in  gun  tubes  is  generally  computed  by 
using  the  classical  Land”  solution  for  thick  vailed  cylinders. 

Thermal  stresses  are  calculated  for  temperature  distributions  that 
are  assumed  to  be  logarithmic  l .  This  approach  does  not  allow  for 
consideration  of  the  effects  of  changes  in  cross  section  nor  for  the 
condition  that  the  internal  pressure  is  applied  over  a  limited  and 
changing  region  of  the  tube.  Neither  does  it  permit  the  considera¬ 
tion  of  transient  temperature  distribution  that  may  occur  in  a 
firing  cycle  [2]',  The  effects  of  non-uniform  conditions  along  the 
axis  and  a  transient  temperature  distribution  changes  the  problem 
from  the  classical  one  dimensional  problem  of  Lame'  to  a  problem  in 
two  dimensions  that  does  not  in  general  have  an  analytical  solution. 

An  approach  with  some  success  for  finding  numerical  solutions  to 
two  dimensional  problems  in  elasticity  is  the  use  of  finite  differ¬ 
ences.  Southwell  [3]  was  the  first  to  do  this  successfully  by  using 
a  relaxation  technique.  More  recently  methods  using  block  relaxation 
and  direct  solution  of  the  system  of  difference  equations  have  been 
used  [4,  5,  6,  7,  8]‘. 

In  this  paper,  a  direct  solution  for  the  system  of  difference 
equations  is  used.  Special  handling  of  the  coefficient  matrix  in 
the  computation  routine  makes  possible  the  handling  of  very  large 
systems  of  equations. 
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w;vBi,0MiKT  or  equations: 

The  differential  equations  governing  tlie  behavior  of  a  gun  tube 
arc  developed  in  terms  of  the  radial  and  axial  displacement  which  are 
designated  u  and  w  respectively. 

The  equations  of  equilibrium  for  circular  symmetry  in  terms  of 
stress  and  with  no  body  forces  arc:  [10] 

3or  +  ®Trz  +  cr“°0  n  (1) 


3icr  z  +  3oz  +  Trz 


Figure  1  shows  a  cross-scction  of  a  typical  tube  with  pressure 
loading.  The  coordinate  axis  and  stress  notations  arc  given. 
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A  and  G  are  Lome  constants  and  a  is  the  coefficient  of  linear 
expansion.  Eacli  is  assumed  to  be  temperature  dependent.  To  and  T 
are  the  initial  and  current  temperatures  respectively. 

The  strains  arc  defined  in  terms  of  displacements  by, 

cr  "  5u»  Eg  *  u,  cz  ■  j)w  +  ez 
Sr  r  Sz 

Yrz  “  1“  +  iii,  e  -  cr  +  Eg  (4) 
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Dimensionless  variables  are  defined  by  the  expressions 

u  *»  a  u 
w  »  a  w 

r  -  H  to  i  C  i  1 

2  B  bn  0  <  n  <  1  (5) 
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where  to  depends  on  Che  bore  radius,  a  is  a  characteristic  radius 
and  b  Che  length  of  the  tube.  Equations  (3),  (4)  and  (3)  are 
substituted  into  equations  (1)  and  (2)4  The  following  two  equations 
In  u  and  w  are  obtained. 
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Equations  (6)  and  (7) ,  must  be  satisfied  throughout  the  tube  and 
are  solved  subject  to  the  appropriate  conditions  on  the  boundary. 


FINITE  DIFFERENCE  ANALOG 

The  finite  difference  analog  for  the  differential  equations  (6) 
and  (7)  is  obtained  by  replacing  derivatives  with  finite  difference 
approximations,  To  minimize  the  truncation  error,  a  central 
difference  approximation  is  used  throughout.  The  dimensionless 
spacing  on  5  is  h  B  1/M,  and  on  n  it  is  k  °  1/M.  The  first, 
second  and  cross  derivatives  at  a  nodal  point  (1,  j)  in  the  open 
region  of  the  tube  cross  section  are: 
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If  approximations  of  the  form  in  equations  (8),  (9)  and  (10)  are 
substituted  into  equations  (6)  and  (7),  two  algebraic  equations  of 
the  following  fora  are  obtained  for  each  nodal  point  within  the 
region. 

Cl,k  "i-l.j-l  +  C2,k  “i-l,J-l  +  c3 ,k  ”i-l,j  +  c4,k  “i-l,J 

+  c5,k  «i-l,j+l  +  c6,k  ni-l,j+l  +  c7,k  ”i,j-l  +  c8,k  Gi,J-l 

+  C9,k  «i,j  +  c10,k  Ci,J  +  cll,k  ”i,j+l  +  c12,k  ui,j+l 
+  C13,k  wi+l, j-1  +  C14,k  ui+l,J~l  +  C15,k  wi+l,j 

+  C16,k  ui+l,j  +  c17,k  ''i+lj+l  +  c18,l;  «i+l,J+l 

•  Cl9,k  k  «  1,  2  (11) 

where  the  coefficients  ^  arc  functions  of  material  properties, 

and  will,  because  of  their  dependence  on  temperature,  be  functions 
of  position.  In  the  solution  of  the  problem  the  material  properties 
are  approximated  by  polynomials  in  temperature  and  their  derivatives 
with  respect  to  the  space  variables  are  obtained  by  using  a  central 
difference  approximation. 

The  set  of  equations  (11)  will  be  equal  in  number  to  twice  the 
number  of  nodal  points  within  the  region.  The  equations  will 
include  values  of  u  and  w  at  each  of  the  nodal  points  inside  the 
region  and  also  values  of  u  and  w  at  nodal  points  on  the  boundary 
line  surrounding  the  region.  At  each  nodal  point  on  the  boundary 
line,  two  algebraic  equations  which  express  the  boundary  conditions 
in  finite  difference  form  arc  written.  When  combined  with  equations 
(11)  for  the  interior  a  solvable  set  of  simultaneous  equations  in 
discrete  values  of  u  and  w  are  obtained. 


BOUNDARY  EQUATION'S  : 

Conditions  at  each  nodal  point  along  the  boundary  will  be  in  the 
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form  of  stress  or  displacement!  and  in  some  instances  it  may  be 
possible  to  define  a  line  of  symmetry  that  can  form  a  boundary  line. 
Boundary  equations  are  generated  for  each  type  of  condition  and  the 
program  is  written  so  that  at  each  boundary  point  any  one  of  the 
three  conditions  can  be  prescribed. 

To  bc6t  accommodate  the  three  possible  conditions  at  the 
boundary  the  nodal  spacing  is  chosen  so  that  the  actual  boundary  of 
the  region  always  lies  midway  between  nodal  lines ,  as  illustrated  by 
figure  2.  Two  boundary  equations  must  fce  written  for  each  of  the 
nodal  points  on  the  outer  perimeter.  The  nature  of  the  equations 
will  be  such  as  to  describe  the  condition  at  the  boundary  that  is 
closest  to  the  nodal  point.  Examples  chosen  from  the  bottom  side  of 

the  region  (shown  in  figure  2)  will  be  illustrated.  Suppose  first 

that  displacements  at  the  nodal  line  i  along  the  bottom  boundary  are 
given  and  are  denoted  by  u.  and  w  .  Then  from  the  central  difference 
approximation,  1 

|  <"l,l  +  ”1.2’  "  ”1  (12«> 

f  tJi.l  +  Sl,2)  ‘“I  (l2b» 


If  the  stresses  at  the  boundary  corresponding  to  the  ic  nodal  line 
are  given,  then  from  equations  (3)  and  (4) 
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The  following  difference  equations  follow, 
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The  symmetry  condition  along  a  boundary  is  given  ty  a  mirror  reflec¬ 
tion  in  that  boundary.  For  the  point  at  the  itn  nodal  line  along 
the  bottom  boundary 

*1,2  "  ”i,l  “  0 

“i,2  +  “i,l  "  0 

Sets  of  equations  similar  to  equations  (12),  (14)  or  (15)  are 
written  for  each  point  along  t he  nodal  boundary  line.  An  exception 
occurs  for  a  corner  point  when  displacement  or  stress  Is  prescribed 
at  the  corner.  For  displacement  given  at  a  corner,  the  boundary 
equation  <*  obtained  by  setting  the  displacement  equal  to  the 
average  of  the  four  neighboring  points.  A  stress  condition  at  a 
corner  is  handled  by  prescribing  the  stiess  on  the  slant  surface  of 
a  triangular  section  with  sides  h/2  by  k/2  as  illustrated  in  figure 
3.  Stresses  on  and  Tng  are  written  in  terms  of  the  boundary 

stresses  which  are  in  turn  approximated  by  finite  differences. 

An  interior  comer  such  as  shown  in  figure  4  is  handled  by 
requiring  the  normal  stress  on  each  boundary  be  satisfied.  The 
shear  stress  in  the  comer  is  ignored.  Each  boundary  equation  is 
of  a  form  that  can  be  written  in  the  general  form  of  equation  (11). 
Thus  all  equations  of  the  region  and  the  boundary  can  be  written 
in  the  format  cf  equation  (11).  Two  equations  arc  written  at  each 
nodal  point,  be  it  interior  or  boundary  point,  and  this  defines  a 
system  of  algebraic  equations  that  is  equal  In  number  to  the 
number  of  discrete  values  of  u  and  w. 


(15  a) 
(15b) 


SOLUTION  OF  ALGEBRAIC  SYSTEM 

The  system  of  algebraic  equations  is  arranged  in  an  order 
prescribed  by  varying  j  and  then  i,  i.c.,  the  first  two  equations 
of  the  set  are  written  for  the  nodal  point  at  (1,  1),  the  next  two 
for  (1,  2)  etc.  Equations  for  points  on  nodal  line  i  will 
in  general  involve  points  on  nodal  lines  at  i  -  1  and  i  +  1.  Thus, 
the  equations  arranged  in  the  above  order  leads  to  a  banded 
coefficient  matrix  on  the  unknowns  u  and  w.  For  executing  the 
solution  of  this  system  on  a  digital  computer  the  matrix  is 
partitioned  on  the  basis  of  values  of  i.  Then,  in  an  elimination 
sequence,  it  is  possible  to  deal  only  with  the  three  consecutive 
sub  matrices  that  correspoivds  to  equations  for  points  on  the  nodal 
lines  i  -  1,  1,  and  i  +  1,  These  matrices  arc  brought  into  core, 
the  elimination  operations  performed  and  then  the  core  updated  for 
the  next  elimination  sequence.  This  enables  a  computer  with 
limited  core  capacity  to  handle  very  large  systems  without  great 
sacrifice  in  time  for  loading  and  unloading  core.  It  also  has 
the  effect  of  limiting  only  the  nodal  spacing  in  the  radial 
direction  and  giving  the  capacity  for  very  large  spacing  in  the 
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axial  or  z  directions,  For  the  IUM  7044  system  with  32k  core 
memory  and  unlimited  disk  storage ,  a  system  with  30  nodal  lines 
in  the  redial  direction,  two  equations  for  each  point,  and  any 
number  in  the  axial  direction  con  be  handled. 

The  computational  sequence  is  handled  in  three  independent 
Btepa,  First  the  coefficient  matrix  is  generated.  This  involves 
computing  the  cofficicnts  for  each  of  the  equations  in  the  set 
and  storing  them  in  a  way  to  take  advantage  of  the  relatively  narrow 
band  width  of  the  matrix.  Next,  the  equations  are  solved  for  the 
unknown  values  of  u  and  w  by  Gaussian  elimination.  The  elimination 
routine  is  carried  out  with  only  a  portion  of  the  coefficient 
matrix  in  core  at  any  time. 

Finally  with  the  displacements  known,  a  separate  program 
evaluates  the  stresses  by  using  central  difference  approximations 
for  the  strains. 


NUMERICAL  EXAMPLE  l 

The  problem  chosen  for  sample  computation  is  a  30MM  gun  tube 
of  very  light  weight  design.  Internal  pressures  are  not  high  but 
the  light  weight  barrel  gets  hot  vary  quickly  in  a  sustained 
burst  and  there  is  real  danger  of  catastrophic  failure.  A  sketch 
of  a  cross-section  of  the  tube  is  given  in  figure  l. 

The  physical  properties  of  the  barrel  are  given  below 
for  ambient  temperature.  These  have  been  taken  as  constants  for 
the  results  reported  here. 

A  ■  1.13  x  107  psi 

C  “  1.73  x  107  psi 

a  -  6.0  x  10"6  in/in/°F 

Computations  have  been  made  for  the  foilwing  cases. 

1.  Tube  at  ambient  temperature  with  internal  pressure  of: 

a.  6000  psi  (projectile  at  21") 

b.  3700  psi  (projectile  at  36") 

2.  Tube  at  8  milliseconds  after  the  20th  round  Is  fired 

with  internal  pressure  of: 

a.  6000  psi 

b.  3700  psi 
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The  transient  temperature  used  for  eases  2  (a)  and  (h)  Is 
given  in  figure  5. 

Results  of  the  stress  computations  arc  given  in  figures  5,  7, 
and  8.  Figure  6  is  a  plot  of  the  three  normal  stresses  versus 
radius  for  cases  (la)  and  (2a)  at  a  spot  21"  from  the  breech. 
Figure  7  shows  results  for  cases  (lb)  and  (lc)  at  29"  from  the 
breech.  Conditions  at  different  positions  along  the  barrel  arc 
similar  in  general  form  but  will  differ  in  magnitude.  This  is 
illustrated  in  figure  8  where  the  tangential  stress  is  plotted 
against  distance  from  the  breech  for  cases  (la)  and  (2a). 

In  general,  the  results  arc  as  expected.  For  the  case  of 
ambient  temperature,  the  values  are  very  nearly  the  sai’ic  as  given 
by  the  classical  theory.  For  the  transient  temperature  condition, 
numerical  values  have  not  been  checked  against  other  methods  of 
calculation  but  the  general  trend  of  the  tangential  stress  seems 
correct.  Perhaps  less  expected  is  the  very  high  axial  stress  in 
the  tube  shown  by  oz  in  the  curves. 
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